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Abstract 

A complete variational treatment is provided for a family of spiked-harmonic 
oscillator Hamiltonians H = — -p* + Bx 2 + 4r {B > 0, A>0), for arbitrary a > 0. 
A compact topological proof is presented that the set S = {ip n } of known exact 
solutions for a = 2 constitutes an orthonormal basis of the Hilbert space £2(0,00). 
Closed-form expressions are derived for the matrix elements of H with respect to 
S . These analytical results, and the inclusion of a further free parameter, facilitate 
optimized variational estimation of the eigenvalues of H to high accuracy. 
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1. Introduction 

A family of quantum Hamiltonians known as spiked harmonic oscillators is given 
by the general Hamiltonian operator 

H = -^+x 2 + —, (1.1) 
dx 2 x a 

acting in the Hilbert space L 2 (0, oo). Eigenfunctions V e £2(0,00) of H satisfy the 
Schrodinger equation 

-ip" + {x 2 + = Eip with V(0) = 0. (1.2) 

The function %l> is an eigenfunction corresponding to the eigenvalue E and the con- 
dition V(0) = is called a Dirichlet boundary condition. The name of the operator 
derives from the graphical shape of the full potential 1 V(x) = x 2 + ^ , which shows 
a pronounced peak near the origin for A > 0. Further, the Hamiltonian Eq.(l.l) is 
characterized by means of two parameters 1 , namely A playing the role of a coupling 
constant and a > determining the degree of the singularity of the potential at the 
origin. Therefore, it has been regarded as a two parameter ( a, A ) problem. Recently 
however, the present authors 2 have studied a more general family of Hamiltonians, 
known as generalized spiked harmonic oscillators, determined by 

H = H Q + ±, H Q = -^- 2 +Bx 2 + 4; (B>0,A>0) (1.3) 

as an operator in L 2 (0, 00) . They argued 2 ~ 5 that the basis, constructed from the exact 
solutions of the singular Goldman and Krivchenkov Hamiltonian H a in Eq.(1.3), 
forms a more effective starting point for a perturbative variational treatment of the 
Hamiltonian (1.1) than the basis of the ordinary harmonic oscillator, as used, for 
example, in earlier works of Aguilera-Navarro et al 6 ~ 7 . The objective of the present 
article is the demonstration of this contention, as well as the establishment of a single 
variational method suitable for the complete set of eigenvalues of the Schrodinger 
equation (1.2). 'Suitable' means that, regardless of the values of the parameters A and 
a , the presented variational treatment remains valid for all discrete energy eigenvalues 
E n , n = 0, 1, 2, . . . and further, is also valid in the TV-dimensional case with an arbitrary 
angular-momentum quantum number I . 
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The article is essentially self-contained, providing verification of the most rele- 
vant results concerning the generalized spiked harmonic oscillators (1.3). The proof 
we here give for the completeness of the set 5" = {ip n } of eigenfunctions of H is 
topological in nature and compact. However, the principal advances of the present 
work over Ref.[2] are twofold. Firstly, we have managed now to derive closed-forms 
for the infinite sums representing the matrix elements of H with respect to S; this, 
in turn, has allowed us to use these matrices to find accurate energy upper bounds 
quickly. Secondly, we have included a free parameter A in the perturbation (1.3) and 
consequently a final minimization of the matrix eigenvalues over A generates even 
more accurate eigenvalue approximations, or allows us to determine the eigenvalues 
to a given accuracy with the use of matrices having smaller dimension D. In Sec. 
2, we give a brief history of the spiked harmonic oscillator problem and indicate its 
relevance. Thereafter we present in Sec. 3 the Gol'dman and Krivchenkov potential 2 
as an exactly solvable model of singular type, and prove that the exact solutions 
of its Schrodinger equation form a complete orthonormal basis of the Hilbert space 
i 2 (0,oo) . This topological fact is accomplished by means of using abscissas of plain 
convergence, absolute convergence as well as of holomorphy of the Laplace Trans- 
form, which abscissas are well discussed in the book of Gustav Doetsch. We further 
develop in Sec. 4 explicit expressions for the matrix elements of the Hamiltonian (1.3) 
for arbitrary values of a subjected to the restriction 7 > § with special attention 
given to a = 2,4 and 6 , for which we have the most complete analytic results. Herein 
7 = 1 + \\f\ + 4-A arises out of the energy expressions for the exact solutions of the 
Gol'dman and Krivchenkov Hamiltonian H , whereas 7 > f follows from the ne- 
cessity of taking inner products of functions of type f(x) = x^^^ a+1 ^ 2 exp(~/3x 2 /2)P(x) 
with one another in the Hilbert space L 2 (0, 00), where P(x) is an arbitrary polyno- 
mial. This is followed by our variational treatment in Sec. 5, wherein we compare in 
detail our method with that of the variational approach of Aguliera-Navarro et al 6 ~ 7 . 
An extension to the N -dimensional case with arbitrary angular momentum number 
/ = 0,1,2,... is provided in Sec. 6. Finally, numerical results with detailed compar- 
isons to pervious work, accompanied by concluding remarks, are given in Sec. 7. The 
difficulty concerning convergence of the eigenvalues approximations calculated from 
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the truncated matrix, in particular the slowness of convergence of the eigenvalues, is 
looked at in depth for the case of a = 4 . It is proved that there is a critical value 
A c = 5/4 such that convergence is slow when A < A c and rapid when A > A c . 

2. Background and brief history 

There are several reasons for the interest in the spiked harmonic oscillator 
Hamiltonian 1-42 . First, it represents the simplest model of certain realistic interac- 
tion potentials in atomic, molecular and nuclear physics, and second, its interesting 
intrinsic properties from the viewpoint of mathematical physics. The potential 
V(x) — x 2 + ^ gives rise to a notably long-ago recognized surprise - i. e. no dominance 
of either of the two interaction potential x 2 and takes place. Thus, for all values 
of A > , the term always adds an infinite repulsive barrier near the origin and 
on the other hand, one can never neglect the x 2 term because its absence demolishes 
the existence of the ground-state energy 43-44 . Consequently 7 , the potential x 2 + ^ 
is like a wide valley extending to oo . Further, the study for small A reveals the 
presence of two different behaviours depending on the value of a . When a < 5/2 , 
the ground state energy has a power series expansion 7 in terms of A. For a > 5/2, 
the regular perturbation theory fails badly for this kind of (supersingular) potential. 
Another interesting observation is that after the perturbation \x~ a is shut off 
(A — > ) , permanent and irreversible (vestigial) effects of the interaction remain 5 . 
The latter effect was first noted by Klauder 8 " 10 who discussed this remarkable 
phenomenon in connection with non-renomalizable covariant quantum field theories. 
In particular, potentials V , which are sufficiently singular cannot be smoothly turned 
off ( A — > ) in the Hamiltonian H = H + XV to restore the free Hamiltonian H . 
This is now know as Klauder's phenomenon 11-13 . Indeed, an obvious approach to 
obtain informations about the spectrum of the Hamiltonian Eq.(l.l) is to regard the 
potential perturbation of the well-known harmonic oscillator Hamiltonian. 

By means of the exact solutions u n of the harmonic oscillator, one attempts to 
investigate the behavior of the spectrum in terms of the unperturbed eigenfunctions 
u n of Hq . Detwiler and Klauder 14 realized that the normal perturbation theory 
already fails badly for the first correction of the ground-state energy (n = 0) for 
a > 3 . Examining the asymptotic behavior of the lowest eigenvalues of H by means 
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of variational arguments, they were able to predict the kind of dependence the 
ground state energy E = E{\) has on the coupling A for A sufficiently small. 

Due to the failure of the Rayleigh-Schrodinger series, a modified perturbation 
theory was required, especially to ascertain the higher-order expressions in the energy 
expansion in terms of powers of A . Harrell 1 modified the Rayleigh-Schrodinger series 
by utilizing the standard (WKB)-approximation technique for the lowest few orders. 
This proving to be quite successful, he continued to developed a special perturbation 
theory, now know as singular perturbation theory, and obtained thereby the first few 
terms of the perturbed A -expansion for different values of a. This turned out to be 
a non-power series expansion and in fact was of exactly the same order as of Detwiler 
and Klauder 14 . Its drawback is however, that it is valid only for sufficiently small 
value of the coupling A - i. e. for A < 1 . 

Since the early work of Detwiler and Klauder and Harrell, spiked harmonic- 
oscillator Hamiltonians Eq.(l.l) have become the subject of intensive study lasting 
over three decades 2-42 . Aguilera-Navarro et al 6 managed to obtain a strong coupling 
perturbative expansion ( A > 2 ) for the ground state energy. Shortly after, Aguilera- 
Navarro and Guardiola 7 attempted to find a path leading from the weak-coupling 
regime due to Harrell 1 to the strong coupling regime 6 . Nonetheless, they failed to 
give a general constructive method to relate the two perturbative regimes for arbi- 
trary values of the exponent a . Special attention was given to certain values of a . 
Aguilera-Navarro et al 17 analysed the Hamiltonian (1.1) for a = 1 around the three 
regions A — » -oo,A — > and A — » oo via Rayleigh-Ritz large-order perturbative ex- 
pansions. Aguilera-Navarro et al 18 studied the particular case of a = 4 by using a 
non-orthonormal set satisfying the correct boundary conditions at x — » oo and x — » . 
Exact and approximate (variational) solutions of the ground state energy of the spiked 
harmonic oscillator problem have also been reported for particular parameter-values 
in the perturbation potential. Aside from these results, special methods have been 
developed to compute the eigenvalues numerically with high precision 30-39 . 

Most of these results however, concern themselves with different approximation 
techniques for the ground-state energy for the problem in one spatial dimension. 
Recently, the N -dimensional case has began to attract the attention of many 
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researchers 3-5,41-42 . In earlier work 2-5 we have pointed out the advantages of basing 
our variational analysis on an exactly soluble model, which has itself a singular 
potential term 2 . We examine a family of generalized spiked harmonic-oscillator 
Hamiltonians Eq.(l.l) in terms of a one-dimensional space variable x (0 < x < oo) 
with eigenfunctions satisfying Dirichlet boundary condition as stated in Eq.(1.3), 
that is to say with wave functions vanishing at the boundaries. The singular 
orthonormal basis, consisting of the set of exact solutions of H , serves as a better 
starting point for the analysis of the Hamiltonian Eq.(l.l). In this paper we use 
these exact solutions of H to provide systematic variational solutions of the spiked 
harmonic oscillator Hamiltonians Eq.(l.l). We first discuss the tools we implement 
herein, namely the Gol'dman and Krivchenkov orthonormal basis. 

3. Orthonormality of confluent hypergeometric series and Singular poten- 
tials 

The purpose of the section is to develop the necessary tools for the variational 
approach we shall present later. Specifically, we derive in this section the orthonormal 
basis most suitable for dealing with the Hamiltonian Eq.(1.3). Although the varia- 
tional method does not per se require the the use of a complete orthonormal basis, 
one can employ in some situations an orthogonal, or even a non-orthogonal 18 , set 
of sufficiently smooth functions to attain very effective results. In such cases, how- 
ever, the variational method shall only have meaning for the particular set of chosen 
functions 18 , and the slightest modification of the original problem shall immediately 
necessitate a new set of functions or at least some modifications. Since we are inter- 
ested in a set of basis functions admissible for the entire parameter range of a and 
A in Eqs.(l.l) and (1.3) for our variational method, we develop our solution by using 
the solution of an exactly solvable singular potential whose singularity coincides with 
that of the spiked harmonic oscillator. Therefore, we use an exactly solvable singular 
Hamiltonian known as the Gol'dman and Krivchenkov Hamiltonian 45-46 



namely tp satisfies the Dirichlet boundary condition. The energy spectrum in terms 
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of the parameters V and a is given by 45 

E n = - a v^jn +1 + 1^/1 + 4F a 2 - 2ayfv^j J, (3.2) 

whereas the exact wavefunctions take the form 

il) n {x) = C n x v e-^ x \F 1 {-n,v+l,^-x 2 ) (3.3) 

with i/ = i(l + Vl + 4Vba 2 ) . In terms of the Pochhammer symbols 

(a) = l, (o) fc = o(o + l)(a + 2) . . . (a + fc - 1) = F( ° + fc) , fc = l,2,... 

expressed in terms of the gamma function V , the previous function 1F1 , known as 
the confluent hypergeometric function 54 , is defined by 

n , \ k 

To simplify the notation, we introduce parameters B = V a~ 2 and A = V a 2 , and 
obtain thereby an exact solution to the one-dimensional Schrodinger equation with 
the singular potential 

A 

V(x) = Bx 2 + —, B > 0, A > 0, (3.5) 
whose energy spectrum is now given in terms of parameters A and B by 

£„ = 2/?(2n + 7 ), n = 0,l,2,..., (3.6) 

wherein j3 = Vb and 7 = 1 + ^Vl + . The wave functions in this case have the form 

il> n (x) = C n x^e-^ x \F 1 {-n;r,Px 2 ) for n = 0,1,2,.... (3.7) 

The constant C n in Eq.(2.7) is determined from the normalization condition 

poo 

/ i>l{x)dx = 1. (3.8) 
Jo 

In order to compute this integral, we need the following lemma, which is a generaliza- 
tion of formula f6 in the appendix of the book Quantum Mechanics by Landau and 
Lifshitz 46 . 

Lemma 1: For 7 > , and m, n = 0, 1,2, . . . 

00 

[ x 2 -<- 1 e-f 3x \F 1 (-n;r,Px 2 )iFi(-m;r,[3x 2 )dx = i^M Sjnn , (3.9) 
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where S mn = for m^n and 1 for m = n. 

Proof: We denote the indefinite integral of the right hand side of Eq.(3.9) by I mn and 
have by means of the series representation Eq.(3.4) of the confluent hypergeometric 
function iFi that 

Further, after resorting to the integral representation of the gamma function 

r oo 

T(x) = / e'H^dt, x>0, 
Jo 

and a change of variables, we obtain for 27 + 2fc + 21 - 1 > that 



In 



_^ n m 

2/97 EE 



(-n) fc (-m)jr(7 + fc + 



^ n r m 

— F F 



(-m)j(7 + fc)j 



(-n) fe r( 7 + fc) 



(7)fc *i 

The finite sum inside the bracket is just the series representation of the terminated 
hypergeometric function 2 F\ , and therefore 

Applying Chu-Vandermonde's theorem 53 on summing the series 2 Fi with unit argu- 
ment, we get 

, _ r( 7 ) " (-fc) m (-„) 



E ^ (3-10) 



wherein we have invoked the identity 



(4)„=(fI if0<n<fc, (3n) 
I if n > k . 

On account of Eq.(3.11), the product {—k) m (—n)k in Eq.(3.10) differs from zero only 
for n = k = m, and thus ( ~ fc)m fc ( ,~" )fc = n\ if n = m or if n ^ m , which terminates 

fe=0 

the proof of our lemma f . 

The reader should realize at this point the connection between the above lemma 
and the ortho normality relation of Laguerre polynomials. In fact, because of the 
relation 54 
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Eq. (3.9) is just another form of the orthonormality condition for Laguerre ployno- 
mials. Thus the normalization constants C n of Eq.(3.7) can now be determined via 
Lemma 1 and the normalization condition Eq.(3.8), and this entails 

c-2 i "!r(7) (312] 

and the normalized wavefunction Eq.(3.7) now reads as follows 



(3.13) 



The alternating coefficients (-1)" were introduced in the definition of ip n (x) to guar- 
antee a smooth transition by means of the identity 

to the odd solutions of the harmonic oscillator problem, namely the case A = in 
Eq.(3.4). From the immediately preceeding we arrive at 

Theorem 1: For ip n (x) defined by Eq.(3.13), the following orthonormality relations 

oo 

J ' ip n (x)ip m {x)dx = 5 mn , TO, 71 = 0,1, 2,... 


hold in L 2 (0, oo) . 

It may seem that the Hamiltonian H is self-adjoint, whence we could assume 
the completeness of the normalized eigenfunctions {ip n (x)}%L from the general the- 
ory of self-adjoint operators. However, we are inhibited from doing this because of 
the domain of definition of H . The Hamiltonian operator H in the Hilbert space 
L 2 (0,oo) cannot have all of L 2 (0,oo) as its domain of definition on account of the pres- 
ence of the second derivative in H . Consequently, the totality C£°(0, oo) of infinitely 
differentiable complex valued functions on (0, oo) with compact support is initially 
assumed to be the domain of definition of H . Cq°(0,oo) lies dense in L 2 (0, oo) and 
H is formally adjoint to itself, but not necessarily self-adjoint. To make H self- 
adjoint, we invoke its Friedrichs extension 47-48 , which is a self-adjoint and exists in 
principle. If we consider H to stand for its Friedrichs extension, then we need to 
know its spectrum a(H ) = a d (H ) u a e (H ) - i.e. the disjoint union of its discrete and 
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essential spectrums. L 2 (0, oo) is the direct sum of the eigenspaces of H if and only 
if, the essential spectrum is empty 49-50 . We establish a e (H ) = by demonstrating 
completeness for the system {ip n (x)}%L in the Hilbert space L 2 (0, oo) of all square 
integrable functions over the interval (0, oo) for H in its Friedrichs extension form, 
by means of the well-known Hahn-Banach theorem. 

Theorem 2: The set of L 2 (0, oo) -functions {^„(a;)}£L defined by Eq.(2.13), is a 
complete orthonormal basis for the Hilbert space L 2 (0,oo) . 

Proof: The orthonormality of {ip n (x)}%L follows from Theorem(l). To prove com- 
pleteness, we proceed as follows. On account of the definition ((3.13)) of ip n (x) in 
terms of the hypergeometric function 

we can express each of the functions n,7, /3x 2 )(n > 0) uniquely as a finite linear 

combination of x 2n (n > 0) as well as conversely. This in terms of the span, which we 
denote by (L.H.) , means that (£.H.)(iFi(-n, 7 ,/3:r 2 )(n > 0)) = (L.H.)(x 2n (n > 0)). We 
multiply each member of these two sets by x 7 ~ 1/2 exp(-/3a; 2 /2) , and thereby obtain a 
linear subspace of the Hilbert space L 2 (0, oo) , whose topological closure satifies 

(L.£f.)(--r-i/2 1 f\(-n, 7 ,/?- 2 )exp(-/3 - 2 /2)(n > 0)) = (L.H.)^- 1 /™" ex p(-/3 - 2 /2)(n > 0)), 

and further, orthogonal complementation yields 

[(L.H.)^ 2 ^ (-n, 7 , P- 2 ) exp(-/3 - 2 /2)(n > 0))]^ = [(L.H.)(-^ 2 + 2n exp(-/3 - 2 /2)(n > 0))]^. 

Now we turn to demonstrating that [( J L. J ff.)(^- 1 / 2 iFi(-n, 7 ,/3- 2 )exp(-/3- 2 /2)(n > 0))] 1 - = 0, 
in other words {C„- 7_1 / 2 iFi(— n,7,/3- 2 )exp(— /^ 2 /2)}£L constitutes an orthonormal basis 
the Hilbert space L 2 (0, oo) with C n 's given by (3.12). Assuming that $<s £2(0,00) with 
$ _L (L.H.)(- 7 - 1 / 2 1 F 1 (-n, 7 ,/3- 2 )exp(-/3- 2 /2)(n> 0)) leads directly to 

p OO 

< .7-V2+2n exp (_ (9 .2 /2 ) I $ >= / ^-1/2+2™ ex p(-/3a ; 2 /2)¥(^) ( ix = (n > 0). 

Jo 

We define 

pOO pOO 

/(«)= / x*- 1 / 2 exv{-sx* /2)&{x)te = l e - st F(t)dt, where F(t) = ^ /2 ~ 3/4 $(\/7)/2. 
Jo Jo 
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Out of the fact that $ e L 2 (0,oo) follows 

/>oo P oo 

/ \e- st F{t)\dt = / l/2F' 2 - 3 ' A \$(Vi)\exp(-{9t8)t)dt 
Jo Jo 




= v/2(5ft s )-7r( 7 ) ||$|| M0iOo) <oo. 

Herewith the abscissa 5 of absolute convergence 51 of the Laplace Integral /(s) = 
£{F} , namely the smallest real number a , such that J °° |e~ s *.F(i)|cft < oo for all 3?s > 5 , 
has the property that & < . In consequence of our substitution x = \fl , we have 
further that 

/•OO 

< •^- 1 / 2+2 ™exp(-/3- 2 /2) | $ >= (-1)" / (-t) n e- {l3/2)t F(t)dt = f (n) ((3/2) = (n > 0). (3.15) 

Jo 

Since /(s) is definitely holomorphic in the half-plane 3?s > , because the abscissa of 
holomorphy \ satisfies x< (3 < a ((3 being the abscissa of ordinary convergence of 
the Laplace Integral 51 C{F}), f(s) shall clearly have a Taylor series expansion 

oo 

/(*) - ]r(n!)-\/«(/3/2)( s - (3/2Y = (\s - (3/2\ < (3/2) (3.16) 

n=0 

in terms of Eq.(3.15) about the point (3/2 > . The radius of convegence of this Taylor 
series is at least a big as (3/2 , since it may happen that the abscissa of holomorphy 51 \ 
of the Laplace Integral f(s) = £{F} is definitely less than the abscissa (3 of ordinary 
convergence, not to mention that of absolute convergence a . Thus we have that 
f(s) = in the open disc \s-/3/2\ < (3/2 of the s complex plane. By means of the identity 
theorem of holomorphic functions 52 , we have that the Laplace Integral f(s) = £{F} = 
in the half-plane of holomorphy 3?s > x ] and more so in the holomorphy domain 
of / , which contains its half-plane of holomorphy. The uniqueness of the Laplace- 
Transform implies that F(t) = l/2f / 2 ~ 3 / 4 $(\/?) = a. e. in t on the interval [0,oo). 
Hence, J °° \^(Vi)\ 2 (2Vt)- 1 dt = / °° \§{x)\ 2 dx = 0, namely $ = a. e. on (0,oo) , which is 
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to say that {C n ■^~ 1 / 2 if\(-n,7;/3- 2 )}$f =0 constitutes a complete orthonormal basis of 
the Hilbert space L 2 (0,oo) f. 

4. The matrix elements for singular potentials 

It is well-known that the Schrodinger equation, even in the one-dimensional case, 
rarely possesses an exact (analytic) solution. Consequently, a multitude of arduous 
numerical techniques have been implemented to ascertain its energy eigenvalues over 
several decades; from among these, we mention matrix diagonalization. Our primary 
aim in this section is to find the matrix elements of the Hamiltonian Eq.(1.3), whose 
calculation is achieved by means of 

Lemma 2: If 27- a > , then for all pairs of non-negative integers and m and n we 
have that 

oc 

(4.1) 

/^- 7 (f)«r( 7 -f) „, a a a 

2 Wn 3F2( ' m ' 7 ~ 2 ' 1 - 2 ; 7 ' 1 - 2 U; 1} - 

Proof: Let I mn denote the infinite integral on the left hand side of Eq.(4.1). Using 
the series representation Eq.(3.4) of the confluent hypergeometric function 1F1 yields 

^^E ^TwT' w [°° x-^+^-'e-^dx. (4.2) 

By resorting to the integral representation of gamma function, we obtain under the 
condition -|+7+fc+/>0 that 



1 (-m)fc(-n); P'i i a 

^ = 2 £S ^WT-^T ^2 + 7 + fc + '> 

(4.3) 



2 ^ ^ ( 7 ) ; /! 



-m)fc 1 



(7)fc M 

This relation we rewrite in terms of the Pochhammer symbols, use r(-f + 7 + £; + /) = 
(-§ + 7 + fc)/ r ( _ § + 7 + fc) , invoke Chu-Vandermonde's theorem and thus write 

W«« "fe (7)«« r( -2 +7 + ^ 

= r(-| + 7 + fc) 2 F 1 (-n, -| + 7 + fc; 7; 1) ( 4 - 4 ) 
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for the finite sum inside the bracket of Eq.(4.3). Consequently, we arrive at 
m ""2 /3 ^; r( 2 +7 + fc) ( 7 )„ ( 7)fc fc! 

_ 1^-7^2) Vr-1 ^ ^ r(f + n-fc)r(-f + 7 + fc) (45) 
2 P (7)n^ J W r( 7 + A)r(|-fc) 

wherein we have used the identity = (-l) fe (™) . We further simplify the expres- 

sion (4.5) by taking note of 

(a- fe) „= ;;- a)fc( t 

(1 - a - n) fc 

in order to justify the relation 



(7)« ^ ( 7 ) fc (l-f -»)**! 



which is valid for all values of a and 7 such that 7 - f > . This completes the proof 
of the lemma f . 

It should be noted that Eq.(4.1) is a generalization of Eq.(3.9), and under the limit 
process of a — > 0+ Eq.(4.1) reduced to Eq.(3.9). In order to prove this, we proceed by 
proving 

OL Ot (X OL 

lim (-)„ 3 F 2 (-m,^- -, 1 - -; 7, 1 - - - n; 1) = n!<5 m „. (4.7) 
a^0+ 2 2 2 2 

This can be demonstrated, by bringing in the series representation of the terminated 
hypergeometric function 3 F 2 , which leads to 

a a a ^ (~m) fc ( 7 - f ) fc (l - f ) fc (-l)T(l - f ) 
( g )„ 3 F 2 (-m, 7 - -, 1 - -; 7 , 1 - - - n; 1) = ^ ______ 

by means of the identities 

a_ r(l-f-n + fc) ( f )„r(l-f-n + fc) 

1 2 Jfe (i- f) _ n r(i-f) 1 J r(i-f) 

Therefore, by taking the limit as a — > 0+ , we can easily see that 

a bm (-)„ 3 F 2( -m, 7--,l-- ;7 ,l---n;l) = g ^TI) 



V 



(-m) fc (-fc) r 
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where we have made use of the relation (1 + fc)_ n = |^)„ • Again, the product 



(-m)fe(-fe)„ leads to the fact that the sum collapses to n\ , for the case 

k=0 

of m = k = n, and otherwise. 

The matrix elements of the Hamiltonian (1.3) are now given by means of Lemma 
(3) in terms of the infinite integral 

oc 

x~° =< ^m\x- a \Tp n >= C m C n J x^- a ~ x e~ 0x \F x (-n, 7, f3x 2 ) 1 F 1 (-m, 7, /3z 2 )^ (4.8) 



Therefore, as consequence of Lemma 2 and Eq.(3.12), the matrix elements now assume 
the explicit forms 



mn 1 j ; ( 7 ) n T( 7 ) V n!m! 3 21 ' 7 2' 2' 7 ' 2' J ' 1 j 

among which the following is of particular interest 



lS (f)»r(7-t)./( 7 ), 



In the case of a being a non- negative even number ( a = 2,4, 6, . . . ), the hyperge- 
ometric function 3 F 2 in Eq.(4.9) may be looked upon as a terminated polynomial of 
degree 1 - § instead of an m -degree polynomial; thus for n > m and a = 2, 4, 6, . . . we 
have that 

3 21 l 2 j ' 7 2' ' 7 ' 2' j ^ S !( 7 ) s (l-f -n) s ' ( ' 

As a result hereof, the matrix elements Eq.(4.9) further simplify into the closed form 
expressions immediately appearing. These are most suitable for computational pur- 
poses as for the case of 7 > and a = , we indeed shall have 

xl n = {\ % n = m (4.12) 

after using Eq.(4.7); as it should have been expected in this specific case. For the case 
of 7 > 1 and a. = 2 , we have from Eq.(4.9) that 



x~ 2 - < 



'(-ir + "/3^V^gf ifn>m, 

\ (4-13) 
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On the other hand, if the case is 7 > 2 and a = 4 , we then have from Eqs.(4.9) and 
(4.11) that 

' (-l) ro+n ^ffei}yife[7(n - m + 1) + 2m] if n > m, 

(4.14) 

(-l)™ + ^«^[7(m - n + 1) + 2n] if m > n. 

As final case that we illustrate, namely 7 > 3 and a = 6 , we point to the fact that 
Eq.(4.9) lets us deduce 



x~ A 



x~ 6 = < 



( I ^m+n0 3 r( 7 -3) / n!( 7 ) m y 
^ ^ 2 T( 7 +2) V m!( 7 )n 

[(2 + n)(l + n) 7 ( 7 + 1) - 2m(l + n)( 7 - 3) (7 + 1) - m(l - m)( 7 - 2) (7 - 3)] if n > m, 



( 1 \m+n 3 r( 7 -3) / m!( 7 )„ 
I ^ 2 T( 7 +2) V nl( 7 ) m X 

I [(2 + m)(l + m)7(7 + 1) - 2n(l + m)( 7 - 3) (7 + 1) - n(l - n)( 7 - 2) (7 - 3)] if m > n. 

(4.15) 

Actually, we can derive similar expressions for all even integers beyond 6 - i. e. 

a = 8, 10, We point to the cases where m > n, where the derivation is achieved 

by interchanging the order of the summation of Eq.(3.2) and applying therefater 
Eq.(4.12) in reversed order. The particular case of A = (or 7 = | ) and B = 1 , that 
is allows us, of course, to recover the result of Aguilera-Navarro et al 6 as a special 
case. 

5. Variational Approach 

In this section we implement the results developed in the previous section to 
calculate the matrix elements of x~ a by means of a complete orthonormal basis. 
Thereby we shall be able to introduce a variational treatment of the spiked harmonic 
oscillator Hamiltonian given by Eq.(l.l). The principle idea is the representation of 
the Hamiltonian Eq.(l.l) as 

<P n 2 A <P r? 2 A , A A ,r ^ 

H =-dxl + BX + -«=-dxl +BX 2 + ^ + ^-^ ^ 

where at a later point, 0) plays the role of an extra degree of freedom, which 
shall be determined through a minimization procedure. Although the idea is not 
complicated, it has many advantages and we mention just a few. First, the range of a 
is no longer restricted and it can be extended as one pleases, provided the condition 
27 > a or more explicitly A > \(a-2) 2 -\ is satisfied. Second, it substantially reduces 
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the number of basis elements required for the computation of the eigenvalues of the 
Hamiltonian Eq.(l.l), even for the intermediate region Awl. Third, it be can be 
adapted effectively and allows easy handling by means of symbolic software such as 
Mathematica. Fourth, the approach of Aguilera-Navarro et al becomes now a special 
case, namely set A = and B = 1 . Fifth, it can be easily extended to the case of N 
dimensions, where the orbital angular momentum number / is arbitrary, with minor 
modifications only, as we shall see in the next section. 

Let i/j(x) be a trial function for Hamiltonian H given by Eq.(5.1), and let us 
suppose that ip(x) is expandable as a finite linear combinations of the basis functions 
i> n {x) as given by Eq.(3.13) - i. e. 

D-l 

^( x ) = ^2 a n tp n (x). (5.2) 

The problem now is to minimize the eigenenergies of Eq.(5.2) with respect to the 
variational parameters a n , n = 0, 1, . . . , D - 1 in the finite dimensional subspace H D 
spanned by the D functions ip Q , ipi, . . . ,tpD-i ■ However, this is equivalent to diagonal- 
izing the Hamiltonian in Eq.(5.2) in the subspace H D . By separating the Hamiltonian 
Eq.(5.2) into two contributions H Q = -d 2 /dx 2 + Bx 2 + Ax~ 2 and Hj = \x~ a - Ax~ 2 , we 
have 

OO 

H mn = J ip m (x)Hip n (x)dx = <ip m \H \ip n > + <ip m \Hi\tp n >, m, n = 0, 1, 2, ...,D- 1. (5.3) 



Since the matrix representation of H is diagonal in the basis {V'nlo , the first term on 
the right-hand side of Eq.(5.3) yields the exact solution of Gol'dman and Krivchenkov 
potential Eq.(3.6) - i. e. 

<iP m \H \Tp n > = 2f3(2n + 7)5 m „ {fi = y/B, 7 = 1 + ^VT+l~A) 

- and the second term on the right-hand side of Eq.(5.3) is given by 

<ip m \H I \tp n > = \<ip m \x~ a \ip n > - A<ip m \x^ 2 \tp n >, 

where <ip m \x~ a \ip n > and <tp m \x^ 2 \tp n > are given by Eqs. (4.9) and (4.13) respectively. 

Two important observations follow from these results. First, the matrix elements 
of the Hamiltonian (5.1) in terms of the Gol'dman and Krivchenkov basis Eq.(3.13) 
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are given explicitly by ( m, n = 0, 1, 2, .., D - 1 , n> m 



y 1+m ^(f)" r( 7 -f),/(7)n(7)r 



F m „ =2/?(2n + 7 )<W + (-l)" +m A/3^ / " v ' V x 

(7)n r(7) V n!m! 



sF 2 (-m, 7 - « 1 - £ ;7 , 1 - „ - J 1) - (-1)—^ ■ l^k 



(5.4) 



diag 


/ -^00 


Hoi 

Hu . 


Hod-i 

H\D-\ 


\ 




\Hd-w 


Hd-ii 


Hd-id- 


J 



2' 2"' 2' ' y ' 7-l\/ m!( 7 )„ 

taken over the D -dimensional subspace spanned by basis (3.13). The expressions in 
Eq.(5.4) are highly suitable for the systematic computer-aided calculations of the 
energy eigenvalues via diagonalization and subsequent minimization of the matrix 



(5.5) 



Secondly, by increasing the matrix dimension D , we can always improve these upper 
energy bounds. In the variational analysis of the ground state energy of the sin- 
gular potential V(x) = x 2 + \x~ a , Aguilera-Navarro et al 6 utilized an orthonormal 
basis of harmonic oscillator eigenfunctions on the interval (0, oo) - i.e. the set of Her- 
mite functions generated by the non-singular harmonic-oscillator potential x 2 . This 
is equivalent to our basis functions for the case B = 1 and A = as mentioned earlier. 
The shortcomings of their approach, however, are as follows. First, validity only holds 
for a < 3 . Second, a huge set of the basis elements is needed to obtain reasonably 
accurate eigenvalues, and this even for the intermediate region (Awl). 

6. The N -Dimensional Case 

In order to extended the scope of our variational analysis to the TV -dimensional 
spiked harmonic oscillator Hamiltonian Eq.(l.l), we first determine the exact solu- 
tions of the TV -dimensional Schodinger equation with a Gol'dman and Krivchenkov 
potential Eq.(3.5). To do this, we notice that the A term of Gol'dman and 
Krivchenkov potential has the dimensions of kinetic energy, such as the term that 
appears in higher-dimensional systems. We therefore may replace A in Eq.(3.6) with 

^^A + A(A + 1), A = Z+i(jV-3), N>2, (6.1) 

and obtain thereby an exact solutions of N -dimensional radial Schrodinger equation 

- & 2 + + Bx2 + = E ^ 1 - (6 ' 2) 
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Such exact solutions are generated from the well-know solutions of harmonic oscillator 
potential by two simple transformations. We first replace the angular momentum I 
in the harmonic oscillator energy expression j3{An + 21 + 3), n = 0,1,2, ... by -\ + 



\J A + (I + i) 2 , and subsequently replace I with A . Thus, the exact eigenvalues of 
JV -dimensional Schrodinger equation with a Gol'dman and Krivchenkov potential are 

£# = 20(2n + 7jv), M = 0,1,2,..., (6.3) 



where (3 = \f~B and -f N = 1 + ^ A + (A + \) 2 , while the exact eigenfunctions are given 
by 

Mx) = (-^"y^^F^"" 6 "*^ 2 ^ 1 ^^^ 5 ^ 2 ^ (6 - 4) 

The matrix elements in this case turn out to be 

<^ m i\x- a \^ nl > = (-l) n+m f3 



n+ mQ % (f )n T(7Af ~ f ) . /(7w)«(7A r )r 



(7)n r( 7A r) V n\m\ ( 6 _5) 
x 3-F2(-m,7jv- |,1- |;7jv,l-n- |;1) N >2. 

Matrix elements for the special cases of a = 2, 4, 6, . . . are obtained by substituting 



in Eqs.(4. 13-15) for 7 the expression 7^ , where 7^ = 1 + y A + (A + ±) 2 . The matrix 
elements of the Hamiltonian Eq.(5.1) now turn out to be very similar to those in 
Eq.(5.4), namely 

H mn = 2t3{2n + lN )5 mn + \<i; ml \x- a \i) nl >-A<il> ml \x- 2 \i; nl >, N>2. (6.6) 

Since the purpose of the present work is to consider the variational analysis of the 
eigenvalues of the Hamiltonian Eq.(l.l) with Dirichlet boundary condition, we re- 
stricted ourself to the case of N > 2 , in order to avoid problems stemming from the 
degeneracy of the spectrum in the case of N = 1 . The one-dimensional case Eq.(5.4) 
on L 2 (0,oo) is recovered by setting D = 3 and I = 0. Further, the recovery of the 
results of Aguilera-Navarro et al is achieved by substituting A = , B = 1 , D = 3 , 
and I = in Eq.(6.6). 

7. Numerical results 



Although the variational method has been used earlier in one form or another, 
the special form it takes in this present article has a specific purpose. First, the 
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demonstration of the validity of an accurate and uniform approximation not only for 
the ground state energy of Eq.(l.l), but also for the entire spectrum in arbitrary 
dimensions with arbitrary angular-momentum number. And second, its applicability 
to general values of the parameters (a, A) of the potential, subjected to the restriction 
a < 27 . The achievement of these purposes is due to the closed forms we have been able 
to obtain for the matrix-elements Eq.(4.9). At this point we note that the variational 
results are upper bounds to the energy levels, in accordance with the variational 
theorem. 

For a < 2 the variational method gives excellent results for arbitrary value of 
the coupling A . The advantage of the minimization process with respect to the pa- 
rameter A is shown in Table (I). Herein we compare the numerical values obtain 
by diagonalizationed for a = 0.5 for different dimensions of the matrix Eq.(5.5) with 
numerical values obtain by the process of diagonalization and minimization with re- 
spect to A . There is interesting observation, which must be noted here. For given 
a < 2 , the dimension of the matrix required for obtaining results of a given accuracy 
depends on the behavior of the coupling A . For A small however, we can get away 
with a matrix Eq.(5.5) of smaller dimension as compared with case where A is large. 
To further illustrate this point, a one-dimension subspace (lxl matrix) is adequate 
for calculating the eigenvalue 3.001 128 (exact to the accuracy quoted) for A = 0.001 
when a = 1 ; but for A = 10 on the other hand, a matrix at least of size 80 x 80 is 
necessary to obtain the eigenvalue 10.577 48 . In the case of a = 1 , we present in Table 
(II) a comparison of our numerical results to some previous work 17 . For a > 2 , the 
situation is completely reversed. 

For a — 2 , the first variational approximation (subspace of dimension one) to the 
ground state eigenvalues of the spiked harmonic oscillator Hamiltonian is 



The minimization of this expression with respect to the parameter A can be easily 
performed and leads to A = A. On substituting this back into Eq.(7.1) we get E = 
2/3(1 + \\/\ + 4A) , which is the exact result quoted in Eq.(3.6) for n = . 

For arbitrary a, the first variational approximation (subspace of dimension 1) of 



E (a = 2) = 2/3 7 + 



A/3 A(3 




(7.1) 



7 — 1 7 — 1 
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the ground state eigenvalues of the Hamiltonian Eq.(l.l) is 

e = mjnflo = mjn{2/3 7 + A/3 t ^ * 1 - -£-}, 7 = l + -vT+4A (7.2) 

When D = 2 , i.e. subspace of dimension 2, the diagonalization can also be performed 
analytically via the secular equation, that is to say by means of the expression 

e± = mjn£ ± = mm{^[(H 00 + H n ) ± ^(#00 - #n) 2 + (2tfoi) 2 ]}. (7.3) 

where e = e_ and ei = e+ . Moreover, one can obtain analytic expressions for upper 
bound eigenvalues 4 in this case. 

Because of the simple formulas for the matrix elements in the cases of a = 4 and 
a = 6 given by (4.14) and (4.15) respectively, the determination of the energy values 
to any desired accuracy has now been reduced to an easy task. A heuristic scheme 
for ascertaining the eigenvalues to any required number of digits is as follows. The 
eigenvalues obtained from successive levels, such as ( 1 x 1, 2 x 2, . . .) , of the truncated 
matrix (5.5) are compared, and the calculation ceases when the successive eigenvalue 
agree with each other up to the prescribed decimal place. It's sufficient, therefore, 
to use the N -dimensional case for the matrix elements Eq.(6.5) for the calculation 
purpose. To recover the 1-dimensional case, we may set D = 3 and I = 0, and to 
cover Aguilera-Navarro et al's results, we set A = 0, B = 1 , D = 3 , and 1 = 0. Table 
(III) illustrate the use of this procedure for the case of a = 4 and A = 1000 in the 
dimensions N = 2 to 10 . 

Another advantage of the variational approach presented herein is the amount of 
information that we get about the spectrum of the Hamiltonian Eq.(l.l) every time we 
compute the eigenvalues via the diagonalization and minimization. Indeed, we obtain 
for an arbitrary matrix (5.5) of size D a set of upper bounds for the eigenvalues 
E Q ,E 1 ,...,E D -i . Each can be improved by either an increase in the dimension of the 
matrix, or by extracting the desired level through the diagonalization and thereafter 
minimizing with respect to parameter A, which is illustrated in Table(IV). 

For a > 2 and small values of A , the variational method is still applicable, however 
the eigenvalues converges very slowly, and an immense number of matrix elements are 
needed in order to obtain accurate results. In principle, the method is still effective 
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as we see from Table (IV), but at the expense of using a huge number of matrix 
elements. On the other hand, we obtained accurate eigenvalues with the use of a 
modest number of matrix elements for A varying between 1000 to 0.01 as Table (V) 
indicates. An interesting explanation of this lies in the following comment. Let us 
look at the perturbed Hamiltonian 

*-£ + *' + £ + 4-?> CM) 

instead of the one given by Eq.(5.1), and solve the matrix eigenvalue problem by 
expanding the determinant in powers of e up to and including e 2 . It is well known 6 
that we shall end up with the perturbation-like formula (after setting e = 1 ) 

E = E a + A<V>o|*- Q |V<o> - A«p \ x -^ > J2 \<M^ a - Ax-^ n >\^ (? 5) 

where the summation on the right-hand side is finite and D is the number of the 
basis functions used. By analyzing this sum, we shall come to understand the slow 
convergence of the eigenvalues calculated by means of the variational approach. For 
the sake of simplicity, let us restrict ourselves to the case of a = 4 . With the aid of 
Eq.(4.13) and Eq.(4.14), the sum in Eq.(7.5) becomes 

L) ■ -4 _ ^-2|„,, ^12 ir2^_ou — (n+1)! ,,^(n + l)! 



\- |<^o|Ax- 4 -Ax- 2 |y> K >| 2 i r 2 ( 7 -2) ^ yti-iy. ^ 

S ^ ^~^H A S^r +2A(A ~ A(7 - 2)) 5 



(n+1)! 



(7.6) 

in this case. For D —> oo , the sums on the right hand side of Eq.(7.6) have closed form 
expressions in terms of hypergeometric functions, particularly 



|<^ |Ax- 4 ~ Ar- 2 |V>»>| 2 = 1 T 2 ( 7 -2) 
„ #0 E n -E 2 7 r 2 ( 7 ) 



E 



A 2 2 f\(3, 1; 7 + 1; 1) + 2A(A - A(j - 2)) 

(7.7) 



x 3^(3, 1, 1; 2, 7 + 1; 1) + (A - A(j - 2)) 2 4 F 3 (3, 1, 1, 1; 2, 2, 7 + 1; 1) 

Herein the conditions of the convergence of hypergeometric functions 53 guarantee 
that the function 2-F1 converges for 7 > 3, the function 3 F 2 converges for 7 > 2 and 
4F3 converges for 7 > 1 . Thus, Eq.(7.7) holds for 7 > 3 or A > 3.75 in general, but for 
a = 4, the matrix elements Eq.(4.9) holds for 7 > 2 or A > 0.75. This demonstrates 
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the difficulties one encounters with the variational method for small value of A . To 
understand this further, we note the relation between A and the parameter A . Thus 
for A small, A is small; on the other hand, A large implies A large. This is evident 
from the first variational approximation, namely Eq.(7.2) yields 

£ o = m in{2 7+I — I A__--A_}, 7 =1 + Ivf+4A (7.8) 

After differentiating Eq.(7.8) with respect to A, it can shown that the relation be- 
tween A and A is implicitly given by the relation 

( 7 - 2)2(4( 7 -!)'-!) 

A M^TV ' ( } 

from which it is quite easy to see that A is an increasing function of 7 for all 7 > 2 . 
Thus for Eq.(7.7) to converge, we must have 7 > 3 - i.e. A > 1.25. This demon- 
strates that for A < 1.25 , the accuracy of the eigenvalue calculations by means of the 
variational method, necessitates a matrix of very large order. 

8. Conclusion 

In this paper we have carried our study of the spiked harmonic-oscillator problem 
further by expressing the Hamiltonian H given by Eq.(l.l) as the perturbation of 
the singular Gol'dman and Krivchenkov Hamiltonian H = + Bx 2 + £ , where 
the expression \x~ a - Ax~ 2 is looked upon as the perturbation term. We have pro- 
vided a compact proof of the fact that the eigenf unctions of order generated by H 
form a suitable singularity-adapted basis for the appropriate Hilbert space of the full 
problem. Our principal results are fourfold: (1) the proof of the completeness of the 
orthogonal set of normalized eigenvalues of H by means of the domains of conveg- 
ence, absolute convergence, holomorphy of the Laplace Transform; (2) the derivation 
of the compact closed form (1.3) for the matrix elements x m " via infinite integrals 
of products of two confluent hypergeometric functions; (3) a variational method ap- 
plied to our pertubation of the Hamiltonian H , which proves to be very effective 
for accurate calculation the eigenvalues of the spiked harmonic oscillator, as the ac- 
companying numerical calculations verify; (4) we show that this variational approach 
is applicable to the entire discrete spectrum for the spiked harmonic oscillator in 
N -dimensions. 
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Table (I) The rate of convergency for the upper bound of the ground state energy 
E of the Schrodinger equation (-A + x 2 + -^)^p = E ip for A = 0.1 and A = 1000. 
E A=a obtained by diagonalization of the D x D matrix elements, E A obtained by 
diagonalization of the D x D matrix elements then minimizing with respect to the 
parameter A . 





A = 0.1 


A = 1000 


Dx D 


e a=q 


E A 


e a=q 


E A 


1 x 1 


3.102 277 


3.102 185 


1025.765 672 


415.934 312 


2x2 


3.102 167 


3.102 149 


746.081 846 


415.932 051 


3x3 


3.102 151 


3.102 143 


642.417 430 


415.890 659 


5x5 


3.102 143 


3.102 141 


549.825 333 


415.889 798 


10 x 10 


3.102 140 


3.102 139 

(Exact) 


461.349 666 


415.889 785 
(Exact) 
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Table (II) A comparison between the results E RP by means of Ricccati-Pade 17 , 
and E R by means of renomalized series 17 and the result of the present work E 
(correct to 7 digit shown) for a = 1 and various values of the coupling A . 



A 


E RP 


E R 


E 


0.001 


3.001 128 


3.001 143 


3.001 128 


0.01 


3.011 276 


3.011 417 


3.011 276 


0.1 


3.110 9 


3.113 386 


3.112 068 


1 


4.057 906 


4.064 649 


4.057 888 


10 


10.577 483 


10.577 825 


10.577 485 



Spiked harmonic oscillators 



page 29 



Table (III) Upper bounds E u for H = -A + x 2 + ^ for dimension N = 2 to 
10, obtained by diagonalization then minimization of the 30 x 30 matrix elements. 
The "exact" values E were obtained by direct numerical integration of Schrodinger's 
equation. 



N 


E u 


E 


2 


21.350 246 


21.350 246 


3 


21.369 463 


21.369 463 


4 


21.427 056 


21.427 056 


5 


21.522 859 


21.522 859 


6 


21.656 596 


21.656 596 


7 


21.827 883 


21.827 883 


8 


22.036 232 


22.036 232 


9 


22.281 057 


22.281 057 


10 


22.561 680 


22.561 680 



Spiked harmonic oscillators 



page 30 



Table (IV) For arbitrary size DxD of the matrix (5.5) for H = -A + x 2 + ±^ , a set 
of upper bounds for the eigenvalues E ,Ei,...,E D -i that can improve by increasing 
the size D . 



1 X 1 


2x2 


3x3 


4x4 


5x5 


6x6 


7x7 


21.427 79 


21.382 12 


21.374 00 


21.370 07 


21.369 72 


21.369 51 


21.369 46 




26.298 42 


26.189 48 


26.166 99 


26.155 44 


26.154 18 


26.153 40 






31.097 17 


30.919 24 


30.878 34 


30.856 56 


30.851 94 








35.834 86 


35.587 50 


35.525 79 


35.492 11 










40.520 33 


40.205 49 


40.121 62 












45.160 79 


44.781 42 














49.762 16 
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Table (V) Upper bounds E u for H = -^ + x 2 + ^- (a = 4 and 6), for small values 
of A , obtained by diagonalization of the DxD matrix elements then minimizing with 
repect to the parameter A . 





A = 0.0025 


A = 0.005 


A = 0.01 


ERef. 


a = 4 


a = 6 


a = 4 


a = 6 


a = 4 


a = 6 


E DK 


3.106 70 


3.353 95 


3.148 39 


3.423 02 


3.205 27 


3.505 74 


E H 


3.037 61 


3.068 22 


3.053 19 


3.081 13 


3.075 22 


3.096 48 


E K 


3.106 81 


3.353 92 


3.148 35 


3.422 88 


3.205 07 


3.505 45 


E p 


3.103 77 


3.343 05 


3.146 64 


3 413 16 


3.204 42 


3.496 88 


E L 


3.106 81 


3.353 92 


3.148 35 


3.422 88 


3.205 07 


3.505 45 


E u 


3.107 95 


3.354 095 


3.149 00 


3.422 95 


3.205 48 


3.505 49 



' K From Detwiler and Klauder (1975). 

H From Harrell (1977). 

K From Killingbeck (1982). Richardson extrapolation. Correct to the 6 digits shown. 

p From Pade approximant technique (W. Solano- Torres et al, 1992). 

L From Lanczos/grid Method (W. Solano- Torres et al, 1992). Correct to the 6 digits 
shown. 

5 From the present work. 
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Table (VI) A comparison between the results E F of Fernandez 24 , and the results 
E HS of Hall et al 29 and the result of the present work e hsk (correct to 7 digit 
shown) for a = 4 and a = 6 and various values of the coupling A . 





a = 4 


a = 6 


A 


EE 


E HS 




E F 




E 


1000 


21.384 46 


21.370 26 


21.369 462 


12.737 60 


12.725 65 


12.718 617 


100 


11.292 41 


11.265 86 


11.265 080 


8.422 60 


8.420 96 


8.413 358 


10 


6.649 78 


6.609 66 


6.606 622 


6.016 40 


6.014 94 


6.003 209 


1 


4.548 79 


4.504 16 


4.494 179 


4.676 88 


4.684 97 


4.659 940 


0.1 


3.626 44 


3.600 44 


3.575 557 


4.019 15 


4.042 84 


3.915 665 


0.01 


3.237 75 


3.249 80 


3.205 486 


3.524 93 


3.580 70 


3.505 492 



